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Abstract. We present an application of a new formalism to treat the quantum 
transport properties of fully interacting nanoscale junctions We consider a model single- 
molecule nanojunction in the presence of two kinds of electron-vibron interactions. In 
terms of the electron density matrix, one interaction is diagonal in the central region 
and the second off-diagonal between the central region and the left electrode. We use 
a non-equilibrium Green's function technique to calculate the system's properties in 
a self-consistent manner. The interaction self-energies are calculated at the Hartree- 
Fock level in the central region and within a dynamical mcan-field-like approach for 
the crossing interaction. Our calculations are performed for different transport regimes 
ranging from the far off-resonance to the quasi-resonant regime, and for a wide range 
of parameters. They show that a non-equilibrium (i.e. bias dependent) dynamical 
(i.e. energy dependent) rcnormalisation is obtained for the contact between the left 
electrode and the central region in the form of a non-equilibrium renormalisation of the 
lead embedding potential. The conductance is affected by the renormalisation of the 
contact: the amplitude of the main resonance peak is modified as well as 'the lineshapc 
of the first vibron side-band. 
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1. Introduction 

The theory of quantum electronic transport in nano-scale devices has evolved rapidly 
over the past decade, as advances in experimental techniques have made it possible 
to measure transport properties down to single-molecule devices. The development 
of accurate theoretical methods for the description of quantum transport at the single- 
molecule level is essential for continued progress in a number of areas including molecular 
electronics, spintronics, and thermoelectrics pQ. 

One of the longstanding problems in quantum charge transport is the establishment 
of a theoretical framework which allows for quantitatively accurate predictions of 
conductance from first principles. The need for methods going beyond the conventional 
approaches, based on equilibrium electronic structure calculations combined with 
Landauer-like elastic scattering, has been clear for a number of years. Only recently 
have more advanced methods to treat electronic interaction appeared (for example see 
Refs. 0EJS]). 

Alternative frameworks to deal with the steady-state or time-dependent transport 
are given by many-body perturbation theory, based on the non-equilibrium (NE) Green's 
function (GF) formalism: in these approaches, the interactions and (initial) correlations 
are taken into account by using conserving approximations for the many-body self- 
energy [a ei ei El [hi uni nn n2i nai m]. 

Other interactions, such as electron-phonon coupling, also play an important role in 
single-molecule quantum transport. These interactions are more important in nanoscale 
systems, as the electronic probability density is concentrated in a small region of space 
and thus normal screening mechanisms are ineffective. Such interactions are also crucial 
in inelastic electron tunneling spectroscopy. Such a technique constitutes an important 
basis for spectroscopy of molecular junctions, yielding insight into the vibrational modes 
(single molecule phonons called vibrons) and ultimately the atomic structure of the 
junction [15]. 

There have been many theoretical investigations focusing on the effect of electron- 
vibron coupling in molecular- and atomic-scale wires [16-51]. In all these studies, the 
interactions have always been considered to be present in the central region (i.e. the 
molecule) only, with the latter connected to two non-interacting terminals. Interactions 
are also assumed not to cross at the contracts between the central region and the leads. 
When electronic interactions are present throughout the system, as within density- 
functional theory calculations, they are treated only at the mean-field level and do not 
allow for any inelastic scattering events. However, such approximations are valid only in 
a very limited number of practical cases. The interactions, in principle, exist throughout 
the entire system. 

In a recent paper we derived a general expression for the current in nano-scale 
junctions with interaction present everywhere in the system [52]. The importance of 
such extended interactions in nano-scale devices has also been addressed, for electron- 
electron interaction, in recently developed approaches such as Refs. [2| 153]. 
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In the present paper, we apply our formalism [52] to a specific model of single- 
molecule nanoj unctions. We focus on a model system in the presence of electron- vibron 
interaction both within the molecule and between the molecule and one of the leads. 
We adopt a quasiparticle-like approach to treat the crossing interactions (i.e. there are 
some restrictions on the components of the self-energy for the crossing interaction) . We 
show how the interaction crossing at one interface of the molecular nanojunctions affects 
the transport properties by renormalising the coupling at the interface in a dynamical 
and bias-dependent manner. 

The paper is organised as follows: In Sec. [2j we briefly recall the main result of our 
expression for the current in fully interacting systems. In Sec. [3j we present the model 
Hamiltonian for the system which includes two kinds of electron-vibron interaction: a 
Holstein-like Hamiltonian combined with a Su-Schrieffer-Heeger-like Hamiltonian. In 
this section, we also describe how the corresponding self-energies are calculated and 
the implications of such approximations on the current expression at the left and right 
interfaces. In Sec. |4j we show how the non-equilibrium dynamical renormalisation 
affects the generalised embedding potential of the left (L) lead, and how in turn this 
affects the conductance of the nanoj unction. We finally conclude and discuss extensions 
of the present work in Sec. |5j 



2. General theory for quantum transport 



We consider a two-terminal device, made of three regions left-central- right, labelled 
L — C — R, in the steady-state regime. In such a device the interaction — which we 
specifically leave undefined (e.g. electron-electron or electron-phonon) — is assumed to 
be well described in terms of the single-particle self-energy £ MB and spreads throughout 
the entire system. 

We use a compact notation for the Green's function G and the self-energy E matrix 
elements M{uS). They are annotated Mc [Ml or Mr) for the elements in the central 
region C (left L, right R region respectively), and Mlc (or Mcl) and Mrc (or Mcr) 
for the elements between region C and region L or R. There are no direct interactions 
between the two electrodes, i.e. E^/ Hi = 0. 

In Refs. [HUE], we showed that for a finite applied bias V the steady-state current 
Il(V) flowing through the left LC interface is given by: 

(lc - r ■-"+Gun' l v+Gm L c -(r i c 
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By definition T,lc(w) = Vlc + ^xc(^) (similarly for the CL components) where Vlc/cl 
are the nominal coupling matrix elements between the L and C regions. g£(w) are 
the GF of the region L renormalised by the interaction inside that region, where 
x = r,a,< stands for the retarded, advanced and lesser GF components respectively. 
For example, for the advanded and retarded components, we have (g r / a (ui))~ 1 = 
(g r / a (u))~ 1 — T^ B ' r ^ a (uj) where all quantities are defined only in the subspace L. 

The second trace in Eq. Q corresponds to inelastic events induced by the 
interaction in the L lead. At equilibrium, because of the detailed balance equation 
E > G < = S < G > , this term vanishes. At non-equilibrium, this is generally not the 
case. However, when a local detailed balance equation holds, i.e. the system is 
locally in a (quasi)equilibrium state, this terms vanishes since one recovers locally 
£ MB,> G < = H ,< G > . Hence, in practice, we do not have to worry about the infinite 
sum in the trace Trm.[...] since there will always be a region/boundary in the left lead 
L beyond which the system is at (quasi)equilibrium. 

The first trace in the current equation Eq. ([!]) corresponds to a generalisation of the 
result of Meir and Wingreen [55]. It encompasses the cases for which the interactions 
are present in the three L, C, R regions as well as in between the L/C and C/R regions. 
It also bears some resemblance to the expression derived by Meir and Wingreen [55] 
when written as: 



jMW = e f dw^ 



(3) 



where we use the standard definitions E c ,<: = — (E^' < )"'' = V CL g^ Vlc — and 
Y, L d a - E§ r = Vcdgl - gl)V L c = iT L . 

By comparing Eq. ([!]) and Eq. (J3J) , we can see that the quantities T LC (Y LC ) and 
T l LC are now playing the role of the L lead self-energy E£ (H r L ) and respectively 
when the interactions cross at the LC interface. 



3. Model for the interaction 



3.1. Hamiltonians 

We consider a single-molecule junction in the presence of electron-vibron interaction 
both inside the central region and crossing at the contacts. We further concentrate 
on a model for the central region which consists of a single molecular level coupled to 
a single vibrational mode. A full description of our methodology, for the interaction 
inside the region C, is provided in Refs. [56l [57] 158] . Moreover, we consider that some 
electron-vibron interaction exists also at one contact (the left L electrode for instance). 
This model typically corresponds to an experiment for a molecule chemisorbed onto a 
surface (the left electrode) with a tunneling barrier to the right R lead. 

In the following model, we consider two kinds of electron-vibron coupling: a local 
coupling in the sense of a Holstein-like coupling of the electron charge density with 
a internal degree of freedom of vibration inside the central region, and an off-diagonal 
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Figure 1. (Color online) Schematic representation of a central scattering region C 
connected to the left L and right R electrodes. Interactions are given by the coupling 
of the region C to the L (R) electrode {V lc /cl an d Vrcicb)i an d by the many- 
body effects within the central region (Sp B ) and at the LC interface (Efc)- E C B 
corresponds to the coupling of an electron with the vibron mode ujq (with coupling 
strength 70), and £jf<§ corresponds to the coupling of a hopping electron with another 
vibron mode u>a (with coupling strength 7^). 



coupling in the sense of a Su-Schrieffer-Heeger-like coupling [5HJ[T7] to another vibration 
mode involving the hopping of an electron between the central C region and the L 
electrode. A schematic representation of the molecular junction is given in Figure [T] 
The Hamiltonian for the region C is 

H c = e d^d + hco a^a + 7o(c^ + a)d^d, (4) 

where d^ (d) creates (annihilates) an electron in the molecular level eo- The electron 
charge density in the molecular level is coupled to the vibration mode of energy uq via 
the coupling constant 70, and a) (a) creates (annihilates) a vibration quantum in the 
vibron mode uq. The central region C is nominally connected to two (left and right) one- 
dimensional tight-binding chains via the hopping integral t L and t R. The corresponding 
electrode a = L,R self-energy is S^(w) = tl a e lka ^ / (3 a with the dispersion relation 
uj = e a + 2(3 a cos(k a (uj)) where e a and (3 a are the tight-binding on-site and off-diagonal 
elements of the electrode chains. 

To describe the electron- vibron interaction existing across the left contact, we 
consider that the hopping integral t 0L is actually dependent on some generalised 
coordinate X. The latter represents either the displacement of the centre-of-mass of the 
molecule or of some chemical group at the end of the molecule link to the L electrode. 
At the lowest order, the matrix element can be linearised as £ol(A) = toL + t' OL X . Hence 
the hopping of an electron from the C region to the L region (and vice versa) is coupled 
to a vibration mode (of energy uja) via the coupling constant 7^4 (itself related to t' 0L ). 
The corresponding Hamiltonian is given by 

H LC = l A (b j + b)(c[d + Sc L ) + u A b% (5) 

where (b) creates (annihilates) a vibration quantum in the vibron mode ua, the 
generalised coordinate is X = y/Kj {2mA^A){b^ + b), and c\ (c^) creates (annihilates) 
an electron in the level El of the L electrode. 
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The Hamiltonians Eq. Q and Eq. ^ are used to calculate the corresponding 
electron self-energies at different orders of the interaction 7 and 7^ using conventional 
non-equilibrium diagrammatic techniques [56, 58J. 

Furthermore, at equilibrium, the whole system has a single and well-defined Fermi 
level /i eq . A finite bias V, applied across the junction, lifts the Fermi levels as 
jj>L,R = Ai eq + Vl,rgV. The fraction of potential drop [60J at the left contact is t/l 
and rj R = i] L — 1 at the right contact, with \i L — fiR = eV and r\ L e [0, 1]. 

3.2. Self-energies for the interactions 

In this paper, we consider different approximations to the treat the interaction inside 
the central region and the interaction crossing at the left contact. First of all, 
calculating exactly the corresponding interaction self-energies is a tremendous task since, 
in principle, they depend on both the phonon progator and the electron Green's functions 
in all the different parts of the system. Hence, for the first application of our formalism, 
we proceed step by step in terms of the increasing complexity of the interaction; and 
we limit ourselves to approximations for the self-energies that are well known and well 
controlled. 

The electron-vibron self-energies in the central region C are calculated within the 
self-consistent Born approximation (i.e. diagrams with one phonon line). The details of 
these calculations have been reported elsewhere [561158] . For the crossing interaction, we 
consider a mean- field-like approach for the electron-vibron coupling at the LC interface. 
Within such an approximation, we can understand the results of the calculations in terms 
of renormalisation of the LC contact. 

Furthermore considering the crossing interaction occuring only at one interface 
permit us to check and test the consistency of our formalism. Indeed, with no interaction 
crossing at the right CR interface, the current Ir is given by the conventional Meir and 
Wingreen formula, i.e. Eq. ^ for the CR interface. In order to have a consistent 
formalism, we need to have current conservation, and such a constraint is best tested 
with crossing interaction at only one interface. The corresponding results are shown in 
detail in |4~T1 

Within mean-field-like approximations, the effects of the crossing interaction 
correspond to a renormalisation of the coupling in a static or a dynamical manner. The 
corresponding self-energies have only retarded and advanced components and 

The extra inelastic effects included in the components E^, L are neglected altogether. 

x 



Hence the Tq quantities defined in Eq. (|2j) become: 

a yvr 

1 C ~ ^CL iJL ^LCi 

^rL.l yir < v-ir 

L C ~ ^CL iJL ^LCi 

with E L ; C = Vlc + Z LC 



L.l y-ir < v- r (^) 
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Correspondingly, the generalised embedding potentials for the left contact defined 
as Y^' x = \E r CL g L Y* r LC \ x [521 [5j] with x = r, a, ^, are now given by 



(7) 



In the static limit, the mean-field approach leads to the Hartree expression for the 
electron- vibron self-energies at the LC interface: 

Ej'> = -2^M, (8) 



Ye'" 


M 


— ^CLdL^LCi 


Ye'" 


M 


_ ^CLiJL^LCi 


Y C L '<( 




= ^cl9l^lc- 



where 



Kc> = -i I ^ G< c {u). (9) 



The self-energy S^' r//a is independent of the energy u and leads to a static 
(nonetheless bias-dependent) renormalisation of the nominal coupling Vcl = Vlc — toL 
between the L and C regions. This NE renormalisation will induce, among other effects, 
a bias-dependent modification of the broadening of the spectral features of the C region. 
We have analysed these effects in details in Ref. [ST]. The renormalisation is such that 
the amplitude of the current is reduced in comparison with the current values obtained 
when the interaction is present only in the central region. The NE static renormalisation 
of the contact is highly non-linear and non-monotonic in function of the applied bias, 
and the larger effects occur at applied biases corresponding to resonance peaks in the 
dynamical conductance. The conductance is also affected by the NE renormalisation 
of the contact, showing asymmetric broadening around the resonance peaks and some 
slight displacement of the peaks at large bias in function of the coupling strengh 7^. 

Beyond the static limit, we can develop a dynamical mean-field-like approach for 
the electron- vibron coupling following Ref. [62]. The retarded self-energy containing 
all orders of the electron-vibron coupling is expressed as a continued fraction as shown 
analytically in Refs. [631 El]- We have already used such an approach to study the 
transport properties of organic molecular wires which are dominated by the propagation 
of polarons [161 DZ] or solitons [18]. The expression for the corresponding self-energy 

= — 2 

Gfc^-wx)" 1 ~ — 2 (10) 

371 



Gl c (u - 2u A ) 



-1 



Gr Lc^ - 3uj a) 



-1 



where G r LC is the retarded component of the off-diagonal GF between the central region 
and the left lead L. Its closed expression is obtained from the corresponding Dyson 
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equation G r LC = [gT,G] r LC : 



G 



LC 



J LC 



Y, 
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J CL 



J LC 



Y, 
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with g r c {<jj) = [to — eq — S 



MB,r/ 
C 



and Y^' r (u) 



Vcr9 r {^)V RC - 



Such a dynamical mean-field-like approach, which corresponds to a quasi-particle 
approach since the self-energy £2cyc*L can ^ e incorporated into a Schrodinger-like 
equation, is expected to affect the transport properties via the non-equilibrium 
dynamical renormalisation of the contacts, i.e. through the generalised embedding 
potentials Y^ ,x {oS) as well as through the corresponding T£(u;) quantities. 

Finally, note that for the lowest-order expansion, the self-energy Y7 LC takes a simple 
form: 



(12) 



4. Results 



We have perfomed calculations, in a self-consistent manner, for many different values 
of the Hamiltonian parameters. We present below the most characteristic results for 
different transport regimes and for different coupling strengths 7a, while the interaction 
in the region C is taken to be in the intermediate coupling regime Uq = 0.2, 7 = 0.15, 
i.e. 7o/w = 0.75. 

The nominal couplings between the central region and the electrodes toL,R, before 
NE renormalisation, are not too large, so that we can discriminate clearly between the 
different vibron side-band peaks in the spectral functions. The values chosen for the 
parameters are typical values for realistic molecular junctions [051 EH]- In the following 
the current is given in units of charge per time, the conductance in units of the quantum 
of conductance G = 2e 2 /H and the bias V and the embedding potential Y c in natural 
units of energy where e = 1 and h = 1. 



4-1- Conserving approximation 

One of the most important physical conditions that our formalism needs to fulfil is the 
constraint of current conservation. We use conserving approximations to calculate the 
interaction self-energies in the central region C and an quasi-particle approximation for 
crossing interaction at the left interface. 



With our choice for Y7 lc (uj) given by Eq. (10), we find that the left lead generalised 

ft ' 

c 



embedding potential Y^' r (u) is renormalised by the crossing interaction. This is clearly 



seen in Figure 2j which shows the imaginary part of Y^ ,r for different transport regimes, 
different coupling at the LC interface and for both equilibrium and non-equilibrium 
conditions. 

At equilibrium, and in the absence of contact renormalisation, the imaginary 

c 



part of Y c ' r is simply the imaginary part of the conventional L lead self-energy 
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Yi r L (uS) = t^ L e %kL<y ^ I 'Pl and corresponds to a semi-elliptic functions with non-zero values 
within the energy range —2(3l < cu < +2(3 L (see dashed lines in figure [2]). The spectral 



weight of 3m Y7 L is given by J du 3m Ti r L (cu) = t% L J du 3m g r L (uj) = — vrtg L . 

In the presence of the renormalised L contact, we obtained a strong deviation from 
the semi-elliptic function as shown by the (red) thin lines and (black) thick lines in 
Figure |2j which correspond respectively to the equilibrium (bias V = 0) and the non- 
equilibrium (bias V = 1.0) conditions. This result indicates a strong reduction of the 
available transport channels in the renormalised L lead embedding potential for regions 
of energy where E^c/lc^) nas non - zero value. 

However, a very important property is that the total spectral weight of 3m Y7 L is 
conserved for all the calculations we have performed, i.e. for all the different transport 
regimes, coupling strength 7^ and ua at the LC interface and all applied bias. We find 
that in all the cases f du 3m Y c ' r (co) = —Ttt\ L (with maximum deviation of 0.05% 
arising from numerical errors in the results shown in Figure [2]). So in this sense, we can 
say that the approximation for E^cyLC ^ s conserving. 

Correspondingly, we have also carefully checked that the current is conserved for all 
the calculations presented in the present paper, i.e. that \I L + I R \/\I L \ ~ \Il + Ir\/\Ir\ < 
10~ 5 , with the current II at the LC interface is given by expression Eq. ([!]) and the 
current Ir at the CR interface does not contain any contact renormalisation, and hence 
is given by a Meir and Wingreen like expression Eq. (|3]). 

4-2. Dynamical non- equilibrium renormalisation 

The dynamical renormalisation of the L lead embedding potential Y^' r presents features 
(dips) at some energies. Qualitatively speaking, there is a form of correlation between 
them and the features that exist in the spectral function of the central region Aq{u)) = 
(GUu)-GUu))/2m. 

Figure [3] shows the spectral function A c (co), rescaled by some factor, for a set 
of parameters corresponding to panel (a) in figure [2j The bottom panel of figure [3] 
shows the imaginary part of Y^' r {oj) from which the semi-elliptic background has been 
subtracted, i.e. AY^' r = Y^' r — Y7 L . Qualitatively speaking, both quantities present 
similar features, shifted in energy by an amount corresponding to uja {oja = 0.2 in the 
calculations) as expected from the definition of Ti T ^^ lc (uj). 

One can then expect a maximum effect of the dynamical renormalisation in an 
integrated quantity such as the current when the features in Y^' r (co) coincide with 
those in the spectral function Ac{uj) within a given bias energy- window. 

We now consider the modification of the dynamical conductance G(V) = dl/dV 
induced by the crossing interaction E^c/lc^) ^ or different transport regimes and for 
different crossing interaction strength. Figures [1] and [5] show G(V) for three different 
transport regimes for weak to strong crossing interaction strength 7^4, and for two values 
of uj A {oj a = 0.2 for Fig. ||and u A = 0.1 Fig. [sj> . 

The dynamical renormalisation of the left contact slightly affects the main resonance 
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Figure 2. (Color online) Imaginary part of the left- lead generalised embedding 
potential Y c ' v (uj). The dashed lines represent the L lead spectral function in the 
absence of renormalisation at the LC interface (i.e. semi-elliptic density of states). 
The thin (red) lines are the spectral functions 3mY c ,' r (w) at equilibrium (V = 0). 
The thick (black) lines are the non-equilibrium spectral functions ^smY c ' (u>) at finite 
bias (V = 1.0). These show a strong reduction of the spectral functions on energy range 
of Aw 1 around £o. Panel (a) Off- resonant regime = 0.5 and strong coupling at 
the interface u>a — 0.20, 70 = 0.28. (b) Off-resonant regime Eq = 0.5 and weak coupling 
at the interface uja — 0.20, 70 = 0.07. (c) Resonant regime £0 = 0.05 and medium 
coupling at the interface lua = 0.20, 70 = 0.14. (d) Intermediate regime £0 = 0.2 and 
medium coupling at the interface uja = 0.10, 70 = 0.07. The other parameters are 
u = 0.20, 7o = 0.15, t QR = t 0L = 0.22, a = 2.0, e a = 0.0, W = 1. 



peak in the conductance, and more importantly the lineshape of the first vibron side- 
band peak above the main conductance peak. The most important effects are obtained 
for the strong crossing interaction strength ja = 0.28. Even if the ratio 'Ja/^a is 
conserved between the calculations shown in Figure [4] and Figure [5j the absolute value 
of 7a is the crucial quantity that governs the effects of the dynamical renormalisation. 
Furthermore, for the different sets of parameters we have considered, our calculations 
show that the lower order expansion for YT LC (see Eq. ()) provides a good approximation 
to the results obtained with a larger number of levels in the continued fraction expansion 




Figure 3. (Color online) Rescaled spectral function proportional to —^stuGq for the 
central region (top panel) and imaginary part of the L lead generalised embedding 
potential relative to the non-renormalised L lead self-energy: AY C ,r = Y c ' r — J7 L 
(bottom panel). The calculations are performed for the same parameters as in panel 
(a) of Figure[2]and for both zero bias and finite bias V = 1. Here //r = and hl = V. 
The qualitative correlations between the features in AYq' t and in the spectral function 
of the C region —^smG r c are clearly observed. 



of T< r LC (see panel (a) in fig. ). 

We now check the effects of varying the nominal coupling t o on the conductance. 
The results are shown in Figure [6] for the off-resonant transport regime with and without 
dynamical renormalisation. With renormalisation of the left contact, the conductance 
values decrease with increasing coupling t 0a to the leads while the conductance peaks 
broaden. This is seen more clearly in the right panel of Figure [6] where the rescale 
conductance G(V)/G ma _ x is plotted such that the main resonance peak has an amplitude 
of 1. 

In the presence of dynamical renormalisation of the L contact, the main 
conductance peak follows the same behaviour as in the absence of renormalisation. 
However, one observes a highly non-linear behaviour of the modification of the first 
vibron side-band peak. The complete and detailed understanding of such modification 
is rather complex, and beyond the scope of the present paper. However it is strongly 
related to the new features in the renormalised L lead embedding potential which depart 
from the non-interacting case. 
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Figure 4. (Color online) Dynamical conductance G(V) for different transport regimes 
and different coupling strength 7a, and oja = 0.20. Panel (a) Resonant regime 
£o = 0.05. Panel (b) Off-resonant regime £o = 0.2. Panel (c) Off-resonant regime 
£o = 0.5. The non-equilibrium dynamical renormalisation of the LC contact affects 
both the main conductance peak and the first vibron side-band peak. The other 
parameters are lj = 0.20, 7 = 0.f5, t 0R = t QL = 0.22, f} a — 2.0, e a = 0.0, i] V = 1. 



5. Conclusion 



We have studied the transport properties through a two-terminal nanoscale device with 
interactions present not only in the central region but also with interaction crossing at 
the interface between the left lead and the central region. To calculate the current 
for such a fully-interacting system, we have used our recently developed quantum 
transport formula [52j based on the NEGF formalism. As a first practical application, 
we have considered a prototypical single-molecule nanojunction with electron-vibron 
interaction. In terms of the electron density matrix, the interaction is diagonal in the 
central region for the first vibron mode and off-diagonal between the central region 
and the left electrode for the second vibron mode. The interaction self-energies are 
calculated in a self-consistent manner using the lowest order Hartree-Fock-like diagram 
in the central region and a quasi-particle (dynamical mean-field-like) approach for the 
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Figure 5. (Color online) Dynamical conductance G(V) for different transport regimes 
and different coupling strength ja, and loa = 0.10. Panel (a) Resonant regime 
£0 = 0.05. Panel (b) Off-resonant regime £0 = 0.2. Panel (c) Off-resonant regime 
£0 = 0.5. The non-equilibrium dynamical renormalisation of the LC contact affects 
both the main conductance peak and the first vibron side-band peak. The other 
parameters are uj — 0.20, 70 = 0.15, £oB = hL = 0.22, (3 a = 2.0, e a = 0.0, rjv = !• 



crossing interaction. Our calculations were performed for different transport regimes 
ranging from the far off-resonance to the quasi-resonant regime, and for a wide range of 
parameter values. 

They show that, for this model, we obtain a non-equilibrium (i.e. bias dependent) 
dynamical (i.e. energy dependent) renormalisation of the generalised embedding 
potential of the L left lead. The renormalisation is such that the amplitude of the 
corresponding spectral function is reduced around the molecular level energy Eq over an 
energy range roughly equal to the energy support of the spectral density of the central 
region C. This corresponds a reduction of the number of transport channels in the left 
lead at given energy, even if the total spectral weight of the L embedding potential is 
conserved. 

The non-equilibrium dynamical renormalisation of the L contact is highly non- 
linear and non-monotonic in function of the applied bias, and the larger effects occur at 
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Figure 6. (Color online) (Left Panel) Dynamical conductance G(V) for the off- 
resonant transport regime (eo = 0.50) and different norminal coupling strength 
ton = £ol to the leads. (Right Panel) Dynamical conductance GiV) rescaled such 
that the amplitude of the main resonance peak is 1 quantum of conductance. The 
other parameters are — 0.2, 7^ = 0.14 (unless otherwise indicated), cjo = 0.20, 
70 = 0.15, /3 a = 2.0, e a = 0.0, r} V = 1. 



applied bias for which features are present in both the spectral function of the central 
region C and the generalised embedding potential of the L lead. The conductance is 
affected by the NE renormalisation of the contact: the amplitude of the main resonance 
peak is modified as well as the lineshape of the first vibron side-band peak. 

Finally, extensions of the present study are now considered to go beyond the quasi- 
particle approach for the crossing interaction self-energy by including as well as the lesser 
and greater components of the LC interface self-energy. This should lead to dynamical 
NE renormalisation of the contact involving inelastic scattering processes. 
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